


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1982-10 


Particle size determination in small solid 
propellant rocket motors using the 
diffractively scattered light method. 


Cramer, Robert Grewelle. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/20296 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist Ser Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ies) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


i ne ' A et 
we Wab ae mC 
_) 8 i 7haF 144 ; rl Sut 
ALL Ree a5 he 1 ae ‘4 Pp ee El 
wih = ara In ae i’ 
“ —_ : vine 7 i i 
‘ (is ~ . tt 
an ‘ J 
A r ail sea tee ; 
7 : A a 7 b et 
A 1A, ea! A i 
7 ; i 4 
af . 
fy thes vee 
“as i TT 
ae : ape tae ro 
uy a [eat = eit ‘ae ih 
ee i ut 
0 ee tie Ara ' 
- a? ’ 1 
Te= ae oh we iy 
’ ta ft a 
so oats | ven 
a7 a. ’ ss 
A i 
' =" 4h 7 i 
4 as Tara i ;" 
ur - spar a 
af r, ‘ uf 4 ie | | 
if J ' 1; 4 
| rit mom 
cae t ft 
seat 7} 
et , 
Ff 
7 pir 
a7 
af ' 
v \ 
to i} / 7 
i 
( oud Roy ' - 
os : 7F \ 
‘e _“. 1 JM 4 i. ' 
ff i 
" a7 a 
af A ‘ , oat 


== 


wae 


a 


ius 


] 1} : 1% 
i aos eI 


LF i | ay) 
| t at if U1 i 
Co ttweae wae 4G 
pant 
5; 


AWN y 






A ve i 


i! 
he 
a a asf 5 - 


maa ite rs ie \ & A 


iy 
ve 






















NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


PARTICLE SIZE DETERMINATION IN SMALL SOLID 
PROGCE REN ~ROCKED MOTORS USING THE 
Dir PRAGitiVEiiy SCATTERED LIGHT METHOD 


by 


Robert Greweile Cramer, Jr. 


Ocvoper 1932 


Nay ysor: D. W. Netzer 





Approved for public release; distribution unlimited 


= 
T206422 





UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) : 


READ INSTRUCTIONS 
REPORT DOCUMENTATION PAGE SEFORE COMPLETING FORM 


S$. TvPE OF REPORT & PERIOO COVERED 
Engineer's Thesis 


October 1982 


- PERFORMING ORG. REPORT NUMBER 




















4. TITLE (and Subtitte) 
Particle Size Determination in Small 
Solid Propellant Rocket Motors Using the 
Diffractiveiy Scattered Light Method 






7. AUTHOR( a) - CONTRACT OR GRANT NUMBER 8) 













F04611-82-x-0008 





Robert Grewelle Cramer, Jr. 









15. PROGRAM ELEMENT PROJECT 


TASK 
AREA & WORK UNIT NUMBERS 





9. PERFORMING ORGANIZATION NAME ANDO ADOORESS 






Naval Postgraduate School 
Monterey, California 93940 


tt CONTROLLING OF FICE NAME ANO ADORESS 
Air Force Rocket Propulsion Laboratory October 1982 | 


MONITORING AGENCY NAME & ADORESE(II different fram Controlling Office) 1B. SECURITY CLASS. (of thie report) 
Unclassified 























ASSItFIC ATION/ DOWNGRADING 


Se, SECL 
EDULE 


$c 





C 
4 


14. OSTRIBUTION STATEMENT (af thle Report) 


Approved for public release; distribution unlimited 


17. DISTRIBUTION STATEMENT (of the ebetcact entered in Black 20, (f different fram Report) 


18. SUPPLEMENTARY NOTES 





19. KEY WOROS (Continue on reveree efde if nececeary and identify by bieck number) 


Solid Propellant Rocket Motors 
Metallized Propellants 
Light Scattering 


20. ABSTRACT (Continue an reveree side tf neceeeary and identify by sieck mamaber) 


A dual beam apparatus was developed which simultaneously 
measured particle size (D337) at the entrance and exit of an 
exhaust nozzle of a small solid propeliant rocket motor. The 
diameters were determined using measurements of diffractively 
scattered laser power spectra. The apparatus was calibrated 
by using spherical glass beads and aluminum oxide powder. 
Measurements were successfully made at both locations. Because 


OD ven 73 1473 Enrtion oF t Nov 481s OBSOLETE 
S/N 0102-014° 6601 | 


UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (Phen Data Entered) 





UNCLASSIFIED 


A ee a EA RE RS ame ee ee 
SecumMmrly CL A6SIFIC ATION OF TIS PAGElWrem Nore Zateved 


of the presence of char agglomerates in the exhaust, con- 
tinued effort is required to improve the grain design in 
order to obtain consistent "across nozzle" data. 


DD Form. 1473 2 


Jan 5 UNCLASSIFIED 
s/N 0102-014-6601 


enim TE s—s—— 
SECUMITY CLASSIFICATION OF THIS PAGE(WRen Data Entered) 


= ee ee ee 


eR cee a A EE ae A 
oe te er 


a DP ~Attn ec T AS TNTCG A EA A te 
ae sa a Se EN On 





Approved for public release; distribution unlimited 


Particle Size Determination in Small Solid 
Propellant Rocket Motors Using the 
Diffractively Scattered Light Method 


by 


ROecememGsewe | le sCramer, Jr. 
Captain, United States Army 
BoS., Washington State University, 1973 


Submitted in partial fulfillment of the 
requirements for the degrees of 


MASEER OF SCEBNCE IN ABRONAUTICAL ENGINEERING 
and 
AERONAUTICAL ENGINEER 
from the 


NAVAL POSTGRADUATE SCHOOL 
October 1982 


Ml 





DUDLEY KNOX LIORARY 


NAVAL POSTGRADUATE SC 
MONTEREY, CALIF, 93949 eee 


ABSTRACT 


A dual beam apparatus was developed which simultaneously 
measured particle size (D,5) at the entrance and exit of an 
exhaust nozzle of a small solid propellant rocket motor. 

The diameters were determined using measurements of diffrac- 
tively scattered laser power spectra. The apparatus was 
calibrated by using spherical glass beads and aluminum oxide 
powder. Measurements were successfully made at both loca- 
tions. Because of the presence of char agglomerates in the 
exhaust, continued effort is required to improve the grain 


weotgd in Order to Obtain consistent “across nozzle" data. 
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i UNTRODUCTAON 


Aluminum is used as a fuel additive in solid propellants 
to increase the specific impulse of the rocket and to sup- 
press high frequency combustion instability. Although there 
is an increase in specific impulse, the specific impulse ef- 
ficiency may actually decrease when compared with the effi- 
ciency of the same base propellant without aluminum [Ref. 1]. 
This performance loss has been largely attributed to incom- 
plete combustion of the aluminum and/or the formation of 
aluminum oxide particles in the motor cavity which leads to 
two-phase flow losses (i.e. velocity and thermal lags between 
the particles and the gas). The magnitude of the losses as- 
sociated with these mechanisms has been reported to be 1% ia 
loss for 10% unburned aluminum and 2-10% loss in = for two- 
phase flow effects [Ref. 2]. These losses are predicted 
analytically and at present there exists no adequate theo- 
retical model relating particle size to propellant and motor 
parameters. Therefore, the accuracy of performance predic- 
tion is very dependent upon the accuracy of particle size 
data obtained experimentally. 

The process by which aluminized fuel additives proceed 
to final combustion products is exceedingly complex. The 
Final particle size is a function of several factors: ori- 


ginal particle size, propellant properties, operating 





environment (pressure, etc.) and the nozzle design and throat 
size [Ref. 2]. 
Kincaid and Derr [Ref. 3] have presented a comprehensive 


description of the process. A summary of their explanation 


is presented here. They divide the rocket motor into five 
zones: (1) at or near the propellant surface (0.5 to 5 mm); 
(2) within the motor cavity; (3) within the converging sec- 
mon of the nozzle; (4) at the nozzle throat; (5) within 
the diverging section of the nozzle. In zone (1), each par- 
ticle is heated as it emerges at the surface. Its metal core 


melts, expands and eventually ruptures its surrounding oxide 
layer. Some of the original metal particles may ignite and 
leave the burning surface immediately. Other particles may 
accumulate on the surface in groups of 100 to greater than 
1000--sometimes a million [Ref. 4]. These groups leave the 
surface either as clumps or burning droplets, called agglom- 
erates. The aluminum oxide particles produced by aluminum 
combustion process have a bimodal size distribution. Fine 
material (1 um) results from vapor phase combustion. This 
Size is generally considered to be invariant and upon com- 
pletion of aluminum combustion amounts to approximately 80% 
by weight of the total condensed phase products. Coarse 
material (100 um) arises from breakup of the residual cap of 
Al,0, on each burning particle. This large cap remains after 
the combustion of the metal has completed. The size of these 


particles is dependent on: propellant composition, pressure, 
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flow environment and the original size of the particle or 
Beohomerate. It is believed that the degree of agglomera- 
tion on the surface will increase for highly aluminized pro- 


pellants. Hence the proportion of coarse particles would 


increase. 

Some oxidation occurs in zone (2}, but the precise mech- 
anism of oxidation is not known. In addition, some growth 
or Al,03 particles by condensation and coagulation occurs in 


this zone. Behavior within zone (2) has been inferred from 
laboratory data. It has been suggested that the behavior in 
zone (2) would be dependent mainly on the size of the parti- 
cles entering the zone. With increased agglomeration in zone 
(1) and short residence time in zone (2) the sareoses courd 
enter the nozzle and be ejected prior to complete combustion. 
Another complication could be settling of the larger particles. 
In zones (3), (4) and (5), size is affected by wall col- 
lisions, growth of particles, breakup of particles, condensa- 
tion and solidification. Zone (3) is a region of accelerating 
gas velocity. When the two-phase flow of zone (2) enters 
this zone, collisions occur which lead to the growth of 
larger particles. As these particles accelerate, their 
break-up is also expected. Because of the change in cross- 
sectional area, the flow must turn and wall collisions occur. 
As the particles enter the throat of the nozzle (zone 4), 
particle breakup has been observed [Ref. 5]. As the flow 


expands through zone (5), it is cooled, and condensation and 
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solidification of the alumina occurs. The problem of two- 
ohase flow losses occurs within zones (3)-(5) when the con- 
densed phase material cannot maintain both velocity and 
thermal equilibrium with the gaseous phase. According to 
theory, condensed material smaller than 1 tm in diameter en- 
tering the nozzle will produce minimal two phase flow losses 
even for highly aluminized solid propellants. 

At this time the Air Force Improved Solid Performance 
Program (SPP) uses a particle size model developed by Hermsen 
(Ref. 6], based on Daz: a weight (or DeBroucker) mean diame- 
ter. This model is based on correlation of a large amount of 
empirical data from many sources, most of which employed par- 
ticle analysis of collected exhausts. Particle break-up and 
fragmentation mechanisms within the nozzle utilize a model 
based on the critical Weber number concept [Ref. 6]. How- 
ever, test cases have indicated that assumptions made in the 
model yield unrealistic particle distributions [Ref. 7]. 

Very little experimental information is available that re- 
lates nozzle geometry and inlet particle size to the behavior 
of the particles as they pass through the nozzle. Light 
scattering techniques were not considered in obtaining par- 
ticle sizes for Hermsen's data base, since they generally 
result in only D35 (volume-surface mean diameter) and are 
considered biased towards small particles. However, this 
technique has several distinct advantages. It is non- 


intrusive and theoretically can be used at any location in 


ee 





or outside of the motor. Particle collection procedures and 
scanning electron microscope (SEM) evaluation can be used to 
determine the particle size distribution and, if this fits a 
log-probability function properly, the measured D3 can be 
directly related to D 43 [Ref. 8]. 

This investigation was part of a larger ongoing effort at 
the Naval Postgraduate School. Its scope was (a) to develop 
and calibrate a diffractively scattered laser power appara- 
tus, (b) to apply the technique at the entrance and exit of- 
the exhaust nozzie of a small rocket motor, and (c) to ob- 
tain "across nozzle" particle size data from various propel- 
lants which can be used to validate or improve upon the 
existing particle size/behavior models used in the Air Force 


Improved Solid Performance Program (SPP). 


is 





Le Dene e hiy hibyY SCATTERED LIGHT METHOD 


A. INTRODUCTION 

The method has its origin in the Mie theory of light 
scattering by spherical particles. Gumprecht and Sliepcevich 
[Refs. 9 and 10] developed from this a theory describing the 
scattering properties of a polydispersion. Using this theory 
@enm, Sliepcevich, and Tribus [Ref. ll] conducted an investi- 
gation into the determination of particle size distributions 
in a polydispersion. This theory was limited to a small op- 
tical depth and there were limits on particle size and re- 
fractive index [Ref. 12]. A more general theory, applicable 
mewcinite optical depths, arbitrarily large particle size and 
arbitrary refractive index, was presented by Dobbins, Crocco 
and Glassman [Ref. 12]. However, their theory was applicable 
oniy to certain types of droplet size distributions. Roberts 
and Webb [Ref. 13] alleviated this problem and demonstrated 
Mae applicability of the "upper-limit distribution function" 
(ULDF) which was introduced by Mugele and Evans [Ref. 14]. 
Using the ULDF Roberts and Webb concluded that the volume- 
surface mean diameter, D5, of a polydispersion could be 
accurately determined from the intensity of diffractively 
scattered light without any knowledge of the general dis- 
tribution type. Dobbins and Jizmagian [Refs. 15 and 16] 


also conciuded that the mean scattering cross section of 
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particles was primarily dependent on the mean diameter Da5, 
and only weakly dependent on the shape of the size distribu- 
tion function when the refractive index is near unity. The 
ratio of scattered light at two forward angles was shown to 
be relatively insensitive to particle refractive index and 
concentration by Hodkinson [Ref. 17]. He also demonstrated 
that this ratio yields a measure of D.,- 
With the above theories as a basis, experimental methods 
have been developed for estimating the volume-surface mean 
diameter, D,., of a polydispersion of particles. Nejad et al 
{[Ref. 18] used this method for measurement of the mean drop- 
let diameter resulting from atomization of a transverse 
liguid jet in a supersonic airstream. Mean diameter, refrac- 
tive index and volume concentration of aerosols have been 
determined by Powell et al [{Ref. 19] using this technique in 
conjunction with light transmission measurements. More re- 
cently the practicality of this method for determination of 
particle sizes in the exhaust of a solid propellant rocket 


has been generally demonstrated by Karagounis [Ref. 20]. 


Bee LECHNIQUE 
The theory for the scattering of light by a single di- 


electric spherical particle of size number, a, gives the fol- 


mowing ©Xpression for the radiant intensity I(8) [{Ref. 12]: 
: ae 
Htc) Mae) 2 | seme) 4m? 
McG u comn) 9 | er | 7 2 ” 
0 (me—2) (mri) 


Jes 





where D is the diameter, @ is the particle size number (7D/\), 
A is the wavelength of incident light, m is the refractive 
index of the scattering media, $8 is the angle measured from 
the forward direction (radians) due to an incident planar 
wave of irradiance Eo and Jy is the Bessel function of first 
fend Of order unity. 

The three terms in the bracket of Equation (1) represent 
the Fraunhofer diffraction, the optical scattering due to 
refraction of the centrally transmitted ray and the optical 
Beoerering due to a grazingly incident ray, in that order. 
Mier eOollowing restrictions apply to Equation (1) [Ref. 12]: 

(1) The incident light must be planar and monochromatic. 
(2) The forward angle @ must be small (98 = sin 89). 


(3) The particle size number (a = 1D/A), and phase shift 
(9 = 2a(m-1)) should be large. 


(4) The distance between particle and observer should be 
large compared to p2/h. 


(5) The particle must be non-absorbing. 

Dobbins et al {Ref. 12] found that for a polydispersion, 
when the expression for intensity of scattering is normalized 
(by dividing by the intensity of diffractively scattered 
iene in the forward direction (98 = 0)), the second and third 
terms in Equation (1) were small and could be ignored. Thus, 
the expression for the normalized integrated intensity of 
forward scattered light, 1(8), due to a polydispersion of 


large particles is given as: 
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| bes ull , 
_1 | Np) p*ap 
ae 


I(6) = — = 2) 
fw (p)p*ap 
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where D, is maximum diameter expected; and N(D) is a distri- 
bution function such that the integral of N(D) over the di- 
ameter interval is the probability of occurrence of particles 
within that interval. Equation (2) is valid when all parti- 
cles are illuminated equally or when the attenuation of the 
incident beam is slight. 

The transmission law for a polydispersion of particles 


immoiven by [Refs. 12, 15, and 21]: 


Be poe lie (- 35 
E, e = exp 5 Q C.,L/D., (3) 
eee by (Ref. 22]: 
HS =) 5 oe (- 20. c t/op (4) 
Ey pee 2 Tae 32 ) 


where Q is the average extinction coefficient; Cc. and Ca are 


eae volume and mass concentration of particles, respectively, 


Mmpeucoldity; L, path length, and 9, particle density. Q and 


D2, are defined as [Ref. 22]: 
© = JQ(D)N(D)D“AD = 
JN (D) D“AD 
y 3 
N 
Dy = L ‘pps (6) 
YN (D)D° AD 
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Dobbins [Ref. 12] recommends that the optical depth (TL) 
be maintained below 1.5 to prevent distortion of the illumina- 
mien profile, 1(9) vs &G. 

Using the ULDF developed by Mugele and Evans [Ref. 14], 
Roberts and Webb [Ref. 13] have concluded that a value of D35 
may be determined from the intensity of diffractively scat- 
tered light from a polydispersion of spherical particles, for 
a wide range of size distributions. They developed a uni- 
versal illumination profile (shown in Figure 5) which is a 
mean curve showing the relation between normalized intensity 
I1(6), and 6, a reduced angle of scattering equal to ™D59/d. 
Using this curve one can determine mean particle size but not 
the distribution of sizes. 

The value of D354 may then be cGetermined as follows. 

(1) The optical apparatus is calibrated by using parti- 
cles of known size as a scattering media, i.e. I(8@) vs 6 
curves are generated and compared to the universal curve. 

(2) The actual particles, in the current investigation 
mie exhaust of the rocket motor, are introaquced and I(6) vs 
6 1s determined. 

(3) Knowing I(8) and 6, 9 is determined from the cali- 
bration curve and/or the universal curve. 


(4) D5 1s determined using the relation: 


ena 


cD] @D| 
=] > 


(7) 
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Ga APPLICABILITY 

Karagounis [Ref. 20] has demonstrated the applicability 
of this technique in estimating the size of Al,0, warticles 
mimehe exhaust of a rocket motor. The restrictive conditions 
for the applicability of equation (2) were satisfied in the 
present experiment in a Similar manner. 

(1) A He-Ne laser with spatial filter and beam expander- 
Sorlimator provided a planar, monochromatic light source. 


(2) The forward angle 9 was less than 6°. 


(3) Large particle size number (>1) a = 7D/i was satis- 
meaceiwicih A = .63208 um, 0.2 um < D < 20 um, then 1<a< 100). 
(4) Large phase shift was provided (>1). (Estimating 


ce |./6 for A102 (ene 2ogedagdel ~ 0 — 100, 1.529 < 152, ) 


(5) The distance between the photodiode array and the 
particles was set at .8 m, which is considerably greater than 
27, = 
mee = 632 um (for a Dre Of 20 “im)-. 
(6) The final condition calling for non-absorbing parti- 
cles is not met exactly, however the absorbing coefficient 


3 


for Al,O, iS approximately 10> times smaller than the scat- 


a3 
tering coefficient for the wavelength, refractive index and 
Size numbers considered here [Ref. 20]. Plass [Refs. 23-25] 

presents a very detailed discussion of this aspect. 
The additional requirement of optical depth, TL < 1.5 is 


satisfied by not applying this technique to situations where 


B/E, < 0.25 (Eq. 4). 
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Iil. EXPERIMENTAL APPARATUS 


A. LIGHT SCATTERING APPARATUS 

A schematic diagram of the final optical set-up is shcwn 
in Figure 1 and photographs are shown in Figure 2. A sche- 
matic of the Data Acquisition System is shown in Figure 3. 
The instruments and optical equipment are listed in Table I. 
A complete description of the Data Acquisition System is pre- 
sented in Reference 26. The optical equipment was mounted 
on two parallel optical benches, one for the beam passing 
through the exhaust and one for the beam passing through the 
motor at the nozzle entrance. Initially a 5 mw Helium-Neon 
laser was used as the light source. This was later replaced 
by a 10 mw Argon laser, for reasons to be discussed later. 
A spatial filter/beam expander located directly in front of 
the light source produced a uniform light beam, one centi- 
meter in diameter. Initially a pinhole with a 2.2 mm diame- 
ter was located in front of the collimating apparatus. This 
pinhole was used to reduce the size of the laser beam in 
order to decrease the measurement volume. Further experi- 
mentation eventually determined that this pinhole was not 
necessary. The intensity of the beam was decreased uSing a 
1 percent neutral density filter to prevent saturation of 
the linear photodiode arrays. The beam was split using a 


50/50 cube beam splitter. One beam continued through the 
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exhaust, the other traveled to a 90 degree prism where it was 
turned parallel to the first beam and directed through the 
motor at the nozzle entrance. Passage of the beam through 
the motor was afforded by mounting circular glass windows in 
the walls of the motor on either side of nozzle entrance 
area. Both beams were passed through identical condensing 
lenses (5 cm in diameter with focal lengths of 500 mm) and 
focused on linear photodiode arrays. The photodiode array 
assemblies each consisted of a linear photodiode array, a 
"mother board" circuit card and an array board. Each as- 
sembly was mounted in a lighttight black box fitted with a 
narrow band pass filter. Each box was mounted on two trans- 
lation stages to provide vertical and horizontal degrees of 
freedom. The linear photodiode array was a self-scanning 
photodiode array consisting of 1024 silicon photodiodes 
mounted in a vertical row on 25 um centers. The aperture 
width was 26 um. The circuits provided an integrated, 
sampled-and~-held output. The “mother board" contained most 
@eeene Circuitry, including the driver/amplifier, clock, 
start, and end of scan functions. The array board contained 
only the circuitry which must be located close to the array 
itself. The output signal was passed through a variable low 
pass filter (for reasons discussed later) before transmission 
to the data acquisition system. The reading of the array was 
accomplished in 36 milliseconds, which corresponds to a rate 


em sO KHZ. 
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Ere oOcks! MOTOR 


The motor was cylindrical, stainless steel with a copper 


nozzle. The chamber was two inches in diameter and two 
inches deep (Figure 4a). The nozzle throat diameter was .245 
inches (for a chamber pressure of 760 psi). An internally 


burning, six-pointed star (Figure 4b) was selected for the 
grain design in an effort to achieve a period of steady state 
pressure in which to take data. The grains were one inch 
long and two inches in diameter with a web thickness of 0.5 
gmenes. Each grain was cut out manually from a one inch 
thick slab of propellant using frabicated shaped cutters. 

The grains were prepared for ignition by applying a thin 
black powder and glue mixture to the inner surface of the 


grain. Electrically ignited BKNO, was discharged from the 


5 
head end of the motor onto the inner surface of the grain, 
thus igniting the black powder, which in turn ignited the 
propellant. A cylindrical stainless steel tube (8" O.D.) was 
mounted at the exhaust of the motor to collect the samples of 
exhaust particles for SEM evaluation. 

Two circular windows were mounted in the walls cf the 
motor on either side of the nozzle entrance area. The win- 
dows were recessed from the motor chamber. To keep the win- 
dows clean, nitrogen was discharged into this recessed area 


and evenly diffused through a sintered metal filter (Figure 


4c). 
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IV. EXPERIMENTAL PROCEDURE 


The apparatus was established in an open test cell of the 
Combustion Laboratory. The output signal of the system con- 
tained a quantity of high frequency noise. After numerous 
unsuccessful attempts to locate and eliminate the source of 
this noise, it was decided to use a variable low pass filter 
to eliminate the noise from the output signal. Subsequent 
examination of the filter effects during the calibration runs 
determined that the filter did not adversely affect the ac- 
curacy of the data readings. After a collimated and filtered 
beam was acquired, the focus was set 1.52 mm above the first 
diode of the array. Care was taken at all times to avoid 
focusing the beam directly on the diodes to prevent damage to 
the array. It was also found necessary to slightly cant the 
focusing lens to prevent interference from reflected light. 

Next, the accuracy of the Diffractively Scattered Light 
Method (DSLM) for obtaining the mean diameter of particles 
was investigated. A comparison was made between experimen- 
tally determined normalized intensity (I(8)) versus non- 
Semensional angle (8) profiles and the “universal" theoretical 
profile (Figure 5). The various profiles were first obtained 
using spherical glass beads of different size ranges. The 
individual profiles for the three size ranges of beads (1-37 


um, 37-44 um and 53-63 um) are shown in Figures 6-8. 
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The profiles for each size range were determined by sus- 
pending the particles in water contained in a Plexiglas box 
(2.5"x3.0" x2.5"). Suspension was maintained using a mag- 
netic stirrer. The first step in this calibration procedure 
was to obtain the zero-scattering profile, that is, the light 
intensity on each diode after passage through the stirred 
water, before particles were introduced. Next the particles 
were introduced and the intensity reading at each diode was 
taken. Examples of intensity versus diode at 1 and 3 kHz 
filter settings are shown in Figures 9 and 10, respectively. 
The analog to digital conversion was accomplished by the data 
acquisition system. (A complete discussion and listing of 
the data acquisition and the reduction procedures discussed 
below are presented separately by Hansen [Ref. 26].) 

After obtaining the two scattering profiles the normal- 
ized intensity I(8) versus angle theta (8) profile was de- 


termined according to the relation: 


_ _78,p 7 79,0_ Dye irons 6,0. 
I(9) = = Fz = 9 FH (8) 
S14 D Gul Cilag cho 
where Tai or Mises are the projected light intensity and re- 
lated voltage at the centerline, and qT, and Va are the light 
iieensity and voltage at an angle, 6 (Figure 11). "o" and "p" 


symbolize no particles and with particles, respectiveiy. The 
First step in determining the normalized intensity was to 


calculate the numerator, AVe, of equation (8) by subtracting 
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the two voltages at corresponding diode numbers. The maximum 
difference (lower limit on diode number) computed was taken 
as the point where the diffractively scattered light inten- 
sity was dominant over the transmitted light intensity. The 
computer determined the diode number where this occurred. An 
upper limit on diode number was taken near the point where 


ee first became zero (Figure 10). A least squares 


Vg,p 9,0 


linear fit was then made to the AV 4 vs diode number data and 


extrapolated to obtain Vol,p = Vane (the forward scattered 
light at 8 = 0). 1(8) vs 8 could then be calculated (equa- 
em@n (8)). 


An example of the relationship of I(8) versus 9 thus pro- 
@mced is shown in Figure 12. At this point the relationship 


of I(6) versus @ is determined using the relationship: 


8 = 9 (9) 





where D3. is the mean diameter of the particle range being 


used in the calibration. 1(8) vs @ can then be plotted (Fig- 
ure 13). An nth order polynomial curve fit was then applied 
to the data over a user specified range. The order of the 


polynomial (also user specified) may be varied from 1 to 10 
depending on the appearance of the data. The range corres- 
ponded to the range of the diodes selected for determination 
of the centerline values as mentioned above. An example of 


the result of this process is shown in Figure 13. The 


oe 





BeSUeeigecurve is then plotted On a semi-log scale for com- 
parison with the universal curve (Figures 5-8). 

The results shown in Figures 5-8 for glass spheres gener- 
ally agree with those of Karagounis [Ref. 20] in that larger 
deviations from the universal curve occur as 8 increases (be- 
yond 6 = 2.5). However, unlike the data of Karagounis, 
agreement with the universal curve was not dependent on 
particle mean diameter. 

Several possible causes for this lack of agreement were 
investigated. First the etfect of filter setting was checked. 
Keeping all other variables constant, data were taken with 
MeeiZ and 3 KHz filter settings. The resultant profiles 
shown in Figure 14 are almost identical. Next the effect of 
Sepa zontal displacement of the array was studied. Again 
keeping all other variables constant, two runs were made. 

The first had the focal point centered on the axis of the 
array. In the second, the focal point was displaced to one 
Side just far enough to minimize non-uniform intensity varia- 
tions. Again the results showed little difference from one 
another (Figure 15). A narrower particle suspension box was 
also used to decrease the depth of the field, but it likewise 
showed little change from the profile using the wider box. 
Signal-to-noise ratio and dark current were two factors which 
could not be completely eliminated and may contribute to the 


lack of agreement at higher 6. 
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Using the same procedure as discussed above, I(98) versus 
68 profiles were prepared for four size ranges of non-spherical 
Al1.,0, powder. The results are shown in Figure 16. The re- 
sults indicate reasonable agreement with the universal curve. 
The experimental profiles had greater slopes than the theo- 
retical profile as they did for the spherical glass beads. 
However, for the Al1,0, the mean size affected the profile. 

The particle size measuring technique was then applied 
eenwactual motor firings. Before each motor firing calibra- 
tion of the recording equipment was conducted. Proper posi- 
tioning of the laser beam on the array was verified. Then 
the zero scattering profile was recorded. The motor was 
fired and the scattering profile was recorded automatically 
at a specified pressure. Chamber vressure was recorded ona 
vVisicorder oscillograph. In the same manner previously de- 
scribed the normalized intensity (1(9)) versus angle (9) was 
determined. D3 values were then calculated over the re- 
corded § range, using both calibration curves and the uni- 


versal curve, from the relation: 


_ Or 
D35 tae (10) 


Fa 


CD 


Residue from the motor exhaust was captured in a stain- 
less steel collection tube. Some of it was dissolved in 
acetone, subjected to an ultrasonic vibrator, allowed to 


settle and then dried in a vacuum oven. This "cleaned" 


a 





sample and an “as collected" sample were examined with an 


SEM. Photographs are shown in Figures 17 and 18. 
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V. RESULTS AND DISCUSSION 


The original goals of this investigation were (a) to 
develop an apparatus to determine D35 from measurements of 
the diffractively scattered laser power spectra using rapid 
scans of iinear photodiode arrays, (b) to apply the apparatus 
at the nozzle exit and through the motor at the exhaust noz- 
zle entrance of a solid propellant rocket motor, and (c) to 
obtain "“across-nozzle" particle data for comparison with the 
existing models for two-phase flow phenomena in exhaust noz- 
zles. The first two goals were met, however, due to experi- 
mental difficulties the last goal was only partially achieved. 

The apparatus developed for this investigation has been 
discussed above. Considerable difficulty was encountered 
during the initial set-up of the computer controlled data 
acquisition and data reducticn procedures. The resulting 
experimental procedures have also been discussed above (Sec- 
tion IV). When the apparatus was applied to actual motor 
firings several additional difficulties occurred. These 
will be discussed below. 

Keeping the windows clear without affecting either motor 
flow or the combustion processes presented the greatest 
challenge. The final solution (Figure 4c) consisted of a 
Sintered metal filter to evenly diffuse the nitrogen into 


the window purge area and a converging nozzle to slightly 
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accelerate the nitrogen flow into the motor. This system 
worked very well at all pressures investigated. 

The light generated during the combustion process within 
the motor (but not at the nozzle exit) was found to produce 
Significant radiation at the wavelength of the He-Ne laser 
mes2c um) initially used as the light source. This pro- 
hibited the measurement of any scattered light from the par- 
ticles. The light source was changed to an argon laser 
(.4880 um) and the problem was alleviated. 

The location of the laser beam in the exhaust jet was 
found to be critical. If directed too far to the rear of 
the nozzle exit, it either missed the jet entirely due to 
the deflection of the jet or the particle number density was 
too low. Both problems resulted in a particle scattering 
profile that was nearly identical to the zero particle scat- 
tering profile. With the current motor the beam was located 
three inches to the rear of the nozzle exit. 

The tests conducted are summarized in Table II. The 
first series of tests was conducted at 775 psia. There 
were two significant results: (1) there was insufficient 
mame through the motor to get a scattering profile; (2) the 
measured D4. in the exhaust was large (30-50 um). The photo- 
graphs in Figures 17 and 18 show the apparent reason for this 


large D to be the presence of a small number of large ir- 


Bz 
regularly shaped agglomerates, which were not Al or Al,0, 


(typically %-5 um spheres). These agglomerates were most 
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Peobasly anhioiteor G6 binder char. Their influence was sig- 
nificant, as even a very few large particles greatly increase 
the measured D,,- The sensitivity of the technique to ir- 
regularly shaped particles has already been demonstrated 
(Figure 16). Repeated tests resulted in different D5, at 
the nozzle exhaust, apparently from different sizes of char 
agglomerates at the time the diode scan was made. 

The tests made at approximately 500 psia (runs 3, 4 and 
5, Table II) were made to determine, (1) if the size of the 
exhaust char agglomerates were sensitive to combustion pres- 
sure, and (2) if the lower chamber pressure would permit 
measurement of scattered light through the motor. [In addi- 
tion, repeated runs were made to acquire the diode data at 
different times during the run to determine whether the data 
was Significantly time dependent. It was found that data 
could be obtained through the motor. Comparisons with the 
universal curve and the small glass bead calibration curve 
(Figure 6) yielded reasonable values for D3. of 16 and 12 um, 
respectively. It was further determined that D25 measured at 
the nozzle exit was sensitive to burn time, as shown by the 
increasing values of D with increased time. And finally, 


Si 


as shown in Figure 17, the spherical A1/A1,0 particles were 


3 


still in the % to 5 um size range. Therefore, the change in 
measured D,, was apparently due to changes in char agglomer- 


ate size during the run. 
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In run 6 an end burning grain was used in an attempt to 
reduce or eliminate the large char agglomerates. The results 
were: (1) a significant reduction in measured D2. (however, 
it was still larger than the A1/A1,0, spheres), and (2) the 
larger optical path with particles (increased optical den- 


sity) precluded the measurement of any scattering data. 


ee 





Viz CONCLUSIONS 


The technique developed for determining the change in 
D3. across a solid propellant rocket motor exhaust nozzle 
appears to work well. Additional effort is required to find 
@a grain configuration and/or an inhibitor that minimizes the 
size of, or eliminates entirely, the char agglomerates. 
This 1S required in order to be able to obtain “across noz- 
Zle" data during a single test. This would allow realization 


Siptne final goal, that is, to obtain data for the validation 


of the two-phase flow models. 
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TABLE I 
EXPERIMENTAL APPARATUS 
Spectra Physics Model 147, Helium-Neon Laser. Power 
5 mW, Random Polarization. 
Power Supply for above, Model 247, 30 watts. 


Laser Beam Expander with spatial filtering, 22 mm aper- 
fice, Oriel, Model 1526. 


Plano-convex lens, 5 cm diameter. 50 cm focal length. (2) 
Neutral density filter, FNG-043, 1%, 2.0 0.D. 

Cube Beam-Splitter, 50/50. 

70> Prism. 


Linear Photodiode Array, EG&G Reticon, G series, 1024 
elements. (2) 


Mother Board Circuit Card, EG&G Reticon, RC1OOB. a) 
Paeay Cil~TCcure Cara, BG&G Reticon, RCIO6. (2) 
Narrow Band Pass Filter. (2) 


DC Power Supply for Photodiode Assembly. 


Variable Low Pass Filter, Krohn-Hite, model 3343. 
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Figure 2; Photographs of Light Scattering Apparatus 


7 





wozshs uot3ytTstnboy ezeq FO oTReWsYoOS 









SAV IS 


LINO OAQIA 24k LVMUY 











l3F AVYUV 
IGOIGOLOHd 


YAILSVW 


SANDING 1a 


YALAdWOOD 


IQOIGOLOHd) 


YOLOW 


ANV114ad0¥d amos aN 


ee 


2-¢ oanbty 


TOYLNOD 


TOYLNOD SYVMaYVH ——— 


SavM1sOS etd — 


YOUMIN 


Y4ILLIIdS WVvae 





Y3SSV1 
3N -2H 


38 





(9zZTS Ten RzoV) 
MeTA-puqg ‘ubtsaq uterzy 


squsuoduojg ACRFOW tp OANhtTyY 


(MOTA OPTS) 
wo3shs ohang CN (9) 





MOPUTM Zaqueuyg 


(q) 





CeO mc 


JOFOW (2) 


TOATUHT | t | 


MOPUTM 


unSvC- 


So 





speog sse[To) Teotaoyds ‘saTtjorg AjtSusqUT. pezt{Tewazon ‘6 SA (6)I ?sg oanbty 


UVa VISHL 


188 *- as ee 18 ° 





BIB" Oa = 
— 4 
rr 
a 
ee 
p—4 
oa 

: ~< 
wn €9-€G -zSpeog SseTHD— --—— \ 
ul PpPRr-LE espeoq sseTH—— -. —— \\\ 
\ 
, wr — >Speog SSeTjS.. ——— ° 
AT Cea peo” ES aI 


DAIND TeSISZAtun 


BT OT 


40 





un £E—-1| “Speeg SSeT[5 [Teotzeuds ‘9 sa (96)I :9 sanbig 


MV VISH! 


188° 


OTA ° 





(ut zt-ez = °*aq) \ 


ul /€-[T ‘-speog sseRTy ss “ 
VAIND TesrAsAtTUy 





BT 


128° 

Ala" = 
- 4 
fvt 
s— 
pod 
—4 
ad 

BI" 

a'r 


41 





wl VP-LE ‘’speod So el.) Teotazayds -9 SA (QE) TI = SANHT yy 


et 


ava VISHL 


100° nn ee eens 





BIB 


(ut ozop = °fa) 
wi prp-L€ espeeq ssetTy—— —— iN 


VAIN [TeSstovATuf) 


B'T 





182 ° 


OTA © 


BT” 


O°T 


ALISNSINI 


42 





ui ¢€9-€6 ‘Speog sseTD TeoTASsyds 49 SA (Q)I *g BaANHTY 


— 


uv ViSHL 


aa 7) OS 


|) ee 


B10" 


(ut pess = °8aq) 
BT” iimGo ec SpoodaSsepo——— ——— i‘ 


Se) wipe ook Uf) 


BT 


120° 

B18" = 
——] 
fvTI 
— 
ly 
of 
—< 

ar” 

aT 


43 





JOYTTA ZH 


T--9potd 


SA OHbeRTOA 


26 eanbt4 





By 


Be * 


By” 


SVL IOA 


ry 
—_ 


44 








ASITTI ZHYA €--Sepotd sa ahe pon 


8 er ee 


-OT S4nbtq 


8 





| po eet 


SOVL IDA 


45 





Azqowosey burza7z34eoSg WybTYT =TT eanbty 


We 


uy 





[6 “s9u] ASIeM JO XSpuT SATRAOeEATOA “Ge°*{T =u 


ail 
@D 


2:297eM AIOY 


{i 


Zt 
x 


Oa sae 108 


Wv38 HaSV] —— pee ‘ r SIX 
PS ee Wo211d0 





Sao iuvd ae I300100LOHd 


46 





Bre” 





oTduexq--9 SA (09)I 


(SNVIGVY) VAHL 


>Z7T eanbty 


B28 ° 


780° 


208 8 


OT 


ALISNSINT CSZPTVAYON 


47 





Ot 


otTduexq--g SA (9)I <€T 


(9) req eqoUL 


—- 


Tq GAINYD 
TetwoudAToOd JapzAO yy" —— —— 


DAINHT Y 


(9) I AQTSusejZUuL paezTTewIo0N 


48 





sbhutqqeg zeqttg Aouenberzyq ZHyY € pue T--@ SA (6)I ‘FI ainbt y 


avd VISHL 
100 "7 je TS A no 8 ie 
A108" \ 


BI” Tower chy €-—- --. —. \ 


JO3TTI ZHY T—- —— \ 
\ 


188 © 


810° 


ALISNAINI 


Bt” 


49 





uostaeduopj AerAyY JO RUDUDDeTdstG TeIUOZTAIOH--g SA (9)I Gt sanbtya 


MV VISHL 


QQ Oo) Q) ™ cO in ~~ 


1) ON aoa 
ph 100° 





a1a° big Ss 
—4 
mM 
<= 
val 
\ a 
— 
X19 oO eee ea ae! 
aL: USD Sepotd —— 
Peas tys 2eeord—~) —— 
Q°T acl 


30 





T3apMOg CO f ty “9 SA (96)I ?9T sanbtg 


— 


ava VISHI 
_ m cD ~ : <t om a “4 S , 
Tae “- Tha 
O18” bin” +S 
es 
71 
—_ 
a 
es ee 7). \ = 
(nee cea ceqmain. ¢9-77——— —— ee a 
SiGe oem cueuinpn= eo = - \ \ \ 
BT” um tg = ““a ‘wi gg-~.—- —— vay a1" 


:zepmog 'o“ty \\ \\ 
Seung Sel ACT) ee \ \\ 


0°T aT B'T 


Sik 








Figure 17: SEM Photograph--Clean Sample (Typical) 
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Figure 18: SEM Photograph--As Collected (Typical) 
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